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ELSEVIER 


Abstract 

The circular duct between the cylinder and displacer serves as a regenerator in free-piston Stirling engines. The cylinder 
wall is fixed and the displacer wall is in reciprocating motion during the steady operation of the engine. The basic equations 
of the working fluid and regenerative duct are derived using the Lagrangian method in terms of the displacement of the 
displacer, so that time does not appear in the equations. A relation is derived between the cylinder and displacer wall tem¬ 
peratures to obtain the initial wall temperature distributions. A computer program is written in fortran and the govern¬ 
ing equations, which include the pressure fluctuations due to the flow reversals, are solved numerically using a finite 
difference method. The results and discussion are presented. © 2002 Elsevier Science Ltd and HR. All rights reserved. 
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Analyse numerique des regenerateurs des moteurs Stirling a 

l’aide de la methode lagrangienne 
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1. Introduction 

In Stirling-cycle machines the working fluid is alter¬ 
natively heated and cooled, and the resulting pressure 
changes in the system perform useful work through 
volume changes. Beale’s free-piston Stirling engine and 
Harwell’s termo-mechanical generator have a duct 
between the cylinder and displacer and it serves as a 
regenerator [1,2]. While the cylinder wall is fixed, the 
displacer wall is in reciprocating motion during the 
operation of the machine. There are many investigations 
on regenerators, but in most cases only their applications 
to the gas turbines or cryogenerators are considered. For 
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the rapidly cycled regenerator of the Stirling-cycle 
machine, pressure fluctuations with time are important. 
The rapidly cycled regenerator theory allows for the 
effects of pressure fluctuations with time by inserting a 
term into the gas energy equation to account for the 
change in internal energy due to pressure pulsation. It is 
essential to know the heat transfer coefficient in order to 
predict the effectiveness of the regenerator. 

In the literature of the theory of rapidly cycled regen¬ 
erators, there are few solutions obtained by approximate 
methods. Assuming a square-wave variation with time 
for the mass flowrate and sawtooth variation with time 
for pressure, Rea and Smith [3] obtained parametric 
solutions for the effectiveness of the regenerator. 
Extending their theory further, Qvale and Smith [4] 
assumed sinusoidal pressure variation and mass flowrate 
with a phase angle between them and obtained an 
approximate solution for the thermal performance of a 
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Nomenclature 

«d 

Velocity of the displacer (m s -1 ) 



u 

Velocity of the gas (ms -1 ) 

A 

Heat transfer area (m 2 ) 

V 

Volume (m 3 ) 

A a 

Cross sectional area of the regenerative duct 

Vs 

Swept volume (m 3 ) 


(nr) 

W 

Amount of work done (J) 

b 

Ratio of heat transfer coefficients (= /j sd r/ d / 

X 

Coordinate of a point on the displacer wall 


Ihy d \) 


from the center of the cylinder (mm) 

c 

Specific heat of the solid (J kg^ 1 K -1 ) 

y 

Displacement of the displacer (mm) 

Cp 

Specific heat at constant pressure (J kg -1 K" 1 ) 

Te 

Entry displacement of displacer for a gas ele¬ 

Cy 

Specific heat at constant volume (J kg -1 K -1 ) 


ment into annular duct (mm) 

d 

Diameter (m) 

Z 

Coordinate of a point on cylinder wall from the 

fc 

Convergence factor 


center of cylinder (mm) 

h 

Heat transfer coefficient (Wm“ 2 K" 1 ) 

q 

Effectiveness of the regenerator 

h s 

Displacement heat transfer coefficient (J nr 3 

(Vq) 

Ineffectiveness of the regenerator 


K- 1 ) 

<p 

Phase angle (rad) 

H 

Enthalpy (J) 

p 

Density (kg m~ 3 ) 

L 

Length of the duct (m) 



m e 

Variation of mass with displacement of dis¬ 

Subscripts 


placer at entry to the duct (kg m -1 ) 

a 

Regenerative duct 

M 

Mass (kg) 

A 

Amplitude 

N 

Velocity ratio ( — ni e RT / A a P) 

C 

Compression volume 

P 

Pressure (Pa) 

d 

Displacer wall 

Q 

Amount of heat transfer (J) 

e 

At entry to the duct 

q s 

Amount of heat transfer per unit displacement 

H 

Expansion volume 


of the displacer (J nr 1 ) 

id 

Ideal 

R 

Gas constant (J kg -1 K -1 ) 

in 

Initial 

S 

Amplitude of displacement movement (mm) 

M 

Mean 

T 

Absolute temperature (K) 

s 

Per unit displacement of the displacer 

T 

- 1 o 

Inlet temperature to expansion volume (K) 

t 

Total 

V 

J o 

Inlet temperature to compression volume (K) 

y 

Cylinder wall 


Stirling-engine regenerator. Kim [5] obtained a closed 
form solution for the amplitude of time-varying instan¬ 
taneous pressure drops by using a sinusoidal pressure 
variation and mass flowrate. 

The most comprehensive approach, known as nodal 
analysis, involves fixed subdivisions of working-fluid 
space in the regenerator. For each subdivision the basic 
equations are solved repeatedly in Eulerian coordinates. 
Using this approach Urieli [6] solved the basic equations 
numerically for each subdivision by a forced-convergence 
technique. Tew et al. [7] analyzed the thermodynamic 
characteristics of a rhombic-drive and a free-piston Stir- 
ling-cycle engine using nodal analysis. Urieli and Kushnir 
[8] showed that this analysis can be utilized in order to 
evaluate the various practical effects of non-ideal regen¬ 
erators, heat exchangers, including heat transfer and 
pressure losses. Organ [9] made a nodal analysis of the 
Stirling-cycle machine and solved the basic equations for 
each subdivision to calculate pressure and gas velocity. 

Rix [10] gave a thermodynamic simulation of a Stir- 
ling-cycle machine by using a Lagrangian formulation. 
In this approach the basic equations are formulated for 


a section of the working fluid by using the gas density, 
velocity and temperature as variables, and cubic spline 
function is used to represent the matrix temperature 
distribution in the regenerator model. 

2. Method of analysis 

In the present analysis using a Lagrangian method, 
the governing equations of the regenerative duct are 
derived in terms of the displacement of the displacer. 
The equations, which include the pressure fluctuations 
due to flow reversals and longitudinal conduction are 
solved numerically by a digital computer using a finite 
difference method. The application of this method to the 
matrix regenerators of the Stirling-cycle machines is 
given in the previous studies [11,12], 

The amount of convective heat transfer to or from a 
gas element based on the displacement of the displacer, 
is defined as 

q s = IhAAT 


( 1 ) 
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where h s is the amount of heat transfer per unit area and 
displacement of the displacer. It is called the displace¬ 
ment heat transfer coefficient in the present analysis and 
is given by 

h s = h/u d (2) 


The ineffectiveness of the regenerator for the com¬ 
pression half-cycle is defined as 

d->/)c= Tu-Tc (3) 

Similarly, for the expansion half-cycle is 


(1 - 0h = 


Th-T q 

Th-Tc 


(4) 


As an approximation the pressure of the system 
assumed to be a linear function of the displacement of 
the displacer and given as 


P=P M +^P A (5) 

Similarly, the masses of the gas in the expansion and 
compression volumes are respectively 


M h = M M h — 4Tah 


2 y 
s 


and 


Mq = Mmc + Mac 


2 y 
S 


( 6 ) 

( 7 ) 


The following additional assumptions are made in 
order to simplify the analysis: 

(a) The displacement heat transfer coefficient is con¬ 
stant. 

(b) The pressure drop across the regenerator is negli¬ 
gible. 

(c) The longitudinal conduction within the cylinder 
and displacer walls and in the gas are negligible. 

(d) The entry temperatures from the expansion and 
compression volumes, 7 h and 7c, are constant. 

(e) There is no radial temperature variation within 
the gas element, displacer or cylinder wall. 

(f) The gas flow in the duct is one-dimensional and 
uniform at any section of the duct. 

(g) The specific heats of the gas, displacer and cylin¬ 
der walls are constant. 

(h) The working gas can be treated as a perfect gas. 
The effect of the power piston is neglected. 

The assumptions of a triangular waveform variation 
in the pressure and the mass of the expansion and com¬ 


pression volumes with the displacement of the displacer 
simplify the analysis of the regenerator more than the 
sinusoidal variations that were assumed by Qvale and 
Smith [4] and Kim [5], and reflect the conditions in a Stir¬ 
ling engine better than the square-wave variation for the 
mass flowrate for pressure with time assumed by Rea and 
Smith [3], Derivation of the basic equations with the 
Lagrangian formulation with respect to the displacement 
of the displacer also simplifies the analysis more than the 
approach used by Rix [10]. Another advantage of this 
approach is that the phase angle between the pressure 
and mass flow does not appear in the equations, and it 
simplifies the analysis of the complicated behavior of the 
machine. The coordinate system used throughout this 
analysis is shown in Fig. 1. 


3. Basic equations 

For the application of the Lagrangian method a gas 
element is considered at the entry to the duct, and 
observations are made on its state as it passes through 
the duct. The gas element is defined by two sections A 
and B, the coordinate of which are functions of two 
independent variables; the instantaneous displacement 
of the displacer, y, and the displacement, y £ at which the 
element under consideration entered the duct, as shown 
in Fig. 2. y and y e are two variables defining the func¬ 
tion z through a relationship of the form 


^ = Aye - y) 


( 8 ) 


At a given displacement, y, the coordinates of the 
points A and B defining the element are 


-a = z(y e ,y) 


( 9 ) 


and 

z B = z(y e + dy e ,y). (10) 

After the displacement, dy the gas element, will 
occupy a different space B'A' in the duct. The coordi¬ 
nates of the points A' and B' will be 

z A ' = z(y e ,y + dy), (11) 


and 


zw = z(y e + dv e , y + dy). (12) 

By using the development of the Taylor series for 
z(y e + dy e , y), neglecting higher order terms and Eq. (9), 
the length of the gas element at the displacement v is 
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3z 1 3 2 z 2 

Z B - r A = Z—fl'Ve + X7T^(«Ve) 
3v e 2 3j2 


(13) 


Owing to the interactions with the surroundings, the 
boundary of the gas element will change, but its mass 
will remain constant. By using a similar method for 
(z B ' — za'X lh e increase in the length of the gas element 
during the displacement, dy of the displacer can be 
written as 


8z S - 


3 2 z 

3v3y e 


dydy e 


(14) 


For the gas element with a constant mass the con¬ 
tinuity equation is written as 


3z m e RT 
3 y e A d P 


(15) 



^y e .y-dy) 


1 

1 


1 


V' 




1 

X 


y 

L. 





(b) 


Fig. 2. Motion of gas element along the regenerator. 


where m e is the mass flowrate of the gas in the duct with 
respect to y at y e . 

To derive the gas energy equation for the gas element 
the first law requires 

8W s = 8Q sy + 8Q sd + dE s (16) 

During the compression half-cycle, the amount of 
heat transferred from the gas element to the cylinder 
and displacer walls is respectively. 

SQ sy = h sy A y (T- T y )dy (17) 


SQ sd = h sd A d (T-T d )dy (18) 

The work done by the gas element per unit displace¬ 
ment of the displacer 

8W s = PAJz (19) 

where Sz is the increase in the length of the gas element 
during the displacement, dy, of the displacer and given by 
Eq. (14). Neglecting the changes in kinetic and potential 
energies, substitution of Eqs. (14), (17), (18) and (19) and 
the value of dE into Eq. (16) will eventually yield 


and 
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A&P + 7tdyh S y(T Ty) 

dy e dy ' 3v e 

+ Jtd d h sd (T- T d )p- + c v m c ^ = 0 
3v e ay 


( 20 ) 


By using Eqs. (5), (7) and (15) with the perfect gas 
law, Eq. (20) will be 


2AaN lt + nd y Nh ^ x + W) ~ T y (z> ~ hT ^ x A 


+ 2 c, 


M a 3T(y ) 
S 3 y 


( 21 ) 


- 0 


3 z A T 
= Ay — 

3.Ve P 

The non-dimensional coefficient Ay in Eq. (28) 
m e RSTc 


A , = 


A&PyyL 


( 28 ) 


(29) 


Similarly, from Eq. (21) the non-dimensional form of 
the gas energy equation can be written as 

ByP k + B 2 (T - T y ) - Bt>(T — T d ) + B 4 % = 0 (30) 

3 y 


where 


b — h sd d d /h S ydy 


and 


N = 


tn e RT 

A- d P 


( 22 ) 


(23) 


From the energy balance for the cylinder wall, the 
cylinder wall energy equation including longitudinal 
conduction, is obtained as 


ndyh S y(T - Ty) + 



3 2 Ty 

3 z 2 


- PyCyA c y 


liTy 

3 y 


(24) 


Similarly the displacer wall energy equation is written 
as 


7r d y h sd (T- T d ) + 



Vd 


3 2 r d 

3 z 2 


— P&CyA 


cy 


3T y 

3 y 


(25) 


The values of non-dimensional coefficients. By, B 2 , B d 


and B 4 are 


II 

(31) 

B _ ndyh sy SNT c 

2P M A a 

(32) 

B _ nd d h, d SNT c 

2T M 4 a 

(33) 

c p M a T c 

4 ~~ PyyAaS 

(34) 


From Eq. (24) the non-dimensional form of the 
cylinder wall energy equation 

/ ~ - \ 3 2 77, dT 

c,(r-r y ) + c 21 J = ¥ (35 ) 

The non-dimensional coefficients Cy and C 2 are 


For simplicity the length of the gas element is 
assumed constant and approximated by 


^. jrdylhyS 

PyCyA a 


Z2 


Ad 

Aa 


dy e 


(26) 


C 2 


ky 

PyV d CyL d 


(36) 

(37) 


4. Nondimensional form of the equations 

The non-dimensional variables are chosen as: 


From Eq. (25) the non-dimensional form of the dis¬ 
placer wall energy equation 



dT 

3 y 


(38) 


z = z/L x = x/L 

y = y/s Te - yjs 

T = T/T c T d = T d /T c , Ty = T y /T c 

P — P/Pm Pa = Pa/Pm 


(27) 


The non-dimensional coefficients Dy and D 2 are 


D { = 


xd d h sd S 

Pd^dAa 


(39) 


By using these non-dimensional variables, the non- 
dimensional form of the continuity equation is 


D 2 


b d 

Pd^AC d L d 


(40) 
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5. Numerical method 

A simplified iterative numerical method is developed 
in order to solve the basic equations of the regenerator. 
This method can be stated in the following steps. 

Step 1. A linear temperature distribution is assumed 
for the displacer wall; 

7d(x) - (T h /Tc - l)x + (Tu/Tc + l)/2 (41) 


Using Eq. (41), Eq. (Al) and (A2) are solved numeri¬ 
cally and the initial temperature distributions of the dis¬ 
placer and cylinder walls, T y (z) and 7d(.v) are obtained. 

By using these initial distributions of T y (z) and Tj(x) 
the temperature distribution of the gas element is 
obtained for a fixed y e during the compression half 
cycle. When the first gas element enters the regenerator 
the initial conditions are; 


y = y e = —1/2, 

(42) 

z= 1/2, 

(43) 

T= Th/Tc 

(44) 

The numerical integration of the gas energy equation 
begins with the above conditions. At the displacement 
of the displacer, y, the conditions for a gas element with 
an entry displacement, y e , will be 

II 

(45) 

- 1/0 (N+l)S^ ^ 
x = 1/2 (y v e ) 

(46) 

<u 

03 1-4 
+ 

<D 

1 

1 

II 

<N 

(47) 


By using the outlet gas temperature from the regen¬ 
erator, the ineffectiveness of the regenerator is calcu¬ 
lated for each half cycle. 

Step 2. The interval of Az, which will ensure the 
numerical stability, is determined using the inequalities 
obtained from Eqs. (35) and (38), and given by 


Cl Ajp e 4 


2C 2 Ay e 
(A zf 


^1 


(48) 


and 


D\ Ay e 4 


2/T Ay c 

(AS ) 2 


^1 


(49) 


By using the gas temperature distribution, the cylin¬ 
der and displacer wall temperature distributions are 


calculated from the cylinder and displacer wall energy 
equations. 

Step 3. For all the A y e intervals selected. Steps 1 and 
2 are repeated for the compression half-cycle. 

Step 4. With the same intervals Ay e and A z, Steps, 1, 
2 and 3 are repeated for the expansion half-cycle. 

Step 5. The net change of the cylinder and displacer 
wall temperatures over the compression and expansion 


half-cycles are written as 


(A7V) C = (tv) c +(4), 

(50) 

( a T)„= (*>)„+(*>), 

(51) 

and 


( 4f '-)c= ( f 0c + (^), 

(52) 

( a T)„= (h) H +<U 

(53) 


The net change of the cylinder and displacer wall 
temperatures over a cycle at each point are expressed as 

A T y = (Ar y ) H +(AT y ) c (54) 

and 

A7d = (ATd)H+(A7d) c (55) 

Initial cylinder and displacer wall temperature dis¬ 
tributions are corrected by these adjustments. 

Step 6. By using the new temperature distribution of 
the cylinder and displacer wall, 

Steps 1-5 are repeated until they converge. The con¬ 
vergence is obtained when the heat transfer over a cycle 
to every point on the cylinder and displacer walls is 
negligible. 

Fourteen points are selected on both the cylinder and 
displacer walls and the amount of heat transfer to and 
from the cylinder and displacer walls are calculated at 
these points for each half cycle. The heat capacity of the 
gas during a blow is much smaller than the heat capacity 
of the walls and the variation of the cylinder and dis¬ 
placer wall temperatures are very small. Therefore, the 
cylinder and displacer wall temperatures are corrected 
after each cycle using Eqs. (54) and (55). 

Although the solution must go through a transient 
before reaching the overall steady-state, the convergence 
is accelerated by a convergence factor to minimize the 
computer time. The net change of the cylinder and dis¬ 
placer wall temperatures over a cycle given by Eqs. (54) 
and (55) are multiplied by this factor. This forced 
convergence technique is similar to that employed by 
Urieli [6] and Rix [10]. The method of the solution is given 
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schematically in Fig. 3 for the compression half-cycle. As 
shown in the figure for each entry displacement, y e , a 
temperature distribution is obtained using the present 
method, and the gas element remained in the duct, is con¬ 
sidered in the analysis. During the flow of gas elements in 
the duct, gas elements are in contact with the same heat 
transfer area on the displacer wall, but with different heat 
transfer area on the cylinder wall as shown in Fig. 3. 

A computer program is written in fortran and by 
using a finite-difference approach the governing equa¬ 
tions are solved on a PC using this numerical method, 
details of which are given by Ataer [13]. 

6. Numerical results 

In the present approach it is also essential to know 
heat transfer coefficients on the cylinder and displacer 
walls, which are assumed constant and calculated from 
the correlations given for periodically reversing flow. 
Air is used as the working fluid. The necessary geome¬ 
trical and thermodynamic data for the numerical solu¬ 
tion of the basic equations are given in Table 1. The 
interval of y e is taken to be 10.0 mm and y to be 1.0 mm. 
The interval of z that will ensure numerical stability is 
calculated using Eqs. (48) and (49). The initial cylinder 
and displacer wall temperature distributions obtained 
from Eqs. (Al) and (A2) are given in Fig. 4. 

For each fixed entry displacement, y e , and with 10.0 mm 
intervals, a gas temperature distribution is obtained, as 
shown in Figs. 5-8. During the compression half cycle 
variation of the gas temperature along the displacer and 
cylinder walls are given in Figs. 5 and 6 respectively for 


different entry displacements. Variation of gas tempera¬ 
ture along the displacer wall, which is in reciprocating 
motion, is shown in Fig. 5. Variation of gas temperature 
along the cylinder wall, which is fixed, is shown in Fig. 6. 
The exit gas temperatures from regenerative duct are 
used to calculate the ineffectiveness of the regenerator. 
The temperature distribution of the gas that remains in 
the regenerative duct during the flow reversal is shown 
by y e = 2 0 mm in Figs. 5 and 6. 

Similarly during the expansion half cycle variation of 
the gas temperature along the displacer and cylinder 
walls are given in Figs. 7 and 8, respectively. The dashed 


Table 1 

Input data for the numerical analysis of the regenerator 


Speed of cycling, co 

400 rpm 

Stroke of the displacer, S 

60 mm 

Compression volume temperature, T c 

387 K 

Temperature difference between the 

150 K 

expansion and compression volumes. 


(Th - T c ) 


Mean velocity of the displacer, u ( \ 

0.71 m/s 

Mean velocity of the air, u 

9.24 m/s 

Velocity ratio, N 

14.29 

Amplitude of the system pressure, T A 

15600 N/m 2 

Amplitude of the mass variation, Mac, M ah 

1.50 g 

Displacement heat transfer coefficient 

383.9 J/m 3 K 

on the cylinder wall, h sy 


Displacement heat transfer coefficient 

414.3 J/m 3 K 

on the displacer wall, h sd 




Fig. 3. Schematic diagram of the method of solution for the compression half-cycle. 
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Fig. 4. Initial cylinder and displacer wall temperature distributions. 



vertical lines in these figures show the location of the gas 
element in the regenerative duct. 

When the steady operation of the engine is reached, 
the temperature distribution of the cylinder and dis¬ 
placer walls are obtained and given in Fig. 9. As seen in 
the figure the temperature distribution of the displacer 
wall slightly deviates from the cylinder wall temperature 
distribution at entry and exit of the regenerative duct. 

As shown in these figures the present analysis computes 
the gas temperature distribution in the regenerative duct 


for each entry displacement of the displacer, and also gives 
cylinder and displacer wall temperatures during steady 
operation of the engine. It simplifies the complicated 
behavior of the free-piston type Stirling engines. 

By using Eqs. (3) and (4), the mean ineffectiveness of 
the regenerator over a cycle is defined as 


(! - ^)m=2 



^oM ^oM 

Tu ~ T c 


(56) 
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For the conditions given in Table 1 and using the data 
given for each fixed entry displacement by Figs. 5-8 the 
mean ineffectiveness of the regenerative duct over a 
cycle is calculated as 16.1%. 

In the present analysis for each fixed entry displace¬ 
ment, y e , the outlet temperature from the regenerative 
duct and ineffectiveness of the regenerator are calculated. 
Examples of the outlet temperatures from the regen¬ 
erative duct are given in Fig. 10, and using the outlet 
temperatures in Eqs. (3) and (4), variation of the inef¬ 
fectiveness with the outlet temperatures for each half¬ 
cycle are also given in Fig. 10. The data obtained from 


the present analysis simplifies the design of the Stirling- 
cycle machines. 

To minimize the computer time the convergence is 
accelerated by a convergence factor. Effect of the con¬ 
vergence factors on the acceleration of the convergence 
is shown in Fig. 11. As shown in the figure the solution 
for f c = 2 converges at cycles higher than 16 cycles, but it 
converges in 10 cycles for f c = 7. An envelope of 
convergence is drawn for the convergence curves as 
shown in Fig. 11. 

The assumption of constant-displacement heat trans¬ 
fer coefficient can be expressed as 
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-120 -too -80 - 60 -40 - 20 0 2 0 40 60 80 100 120 


z mm 

Fig. 8. Gas temperature distribution along the cylinder wall during the compression half-cycle. 



BO DO » 40 » -20 - 50 -80 -110 -00 


z,x mm 

Fig. 9. The cylinder and displacer wall temperature distributions. 



(57) 


and the relation between the gas and the displacer velo¬ 
city approximated as 

Aa 

ii = —— u L \ (58) 


The velocity of the gas appears with the power of 
—0.2 in the expression of displacement heat transfer 
coefficient and the error, caused by the assumption, is 
smaller than in the other theories based on the constant 
heat transfer coefficient 

In the analysis the entry temperatures from the 
expansion and compression volumes, T H and Tc, are 
assumed to be constant. Neglecting the effect of the 
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Fig. 10. Variation of ineffectiveness and outlet temperatures with the entry displacement of the displacer. 



Fig. 11. Effect of different factors on the acceleration of the convergence. 


piston, the volume variation of the expansion space with 
the displacement of the displacer is 


Vh — Fdh + 



(59) 


where F DH and V s are the dead and swept volumes of 
the engine, respectively. Substituting Eqs. (5) and (6) 


and (59) into the perfect gas law, and expanding the 
term [1 — (2y/S)M\y[/Mym] as a binomial, it can be 
approximated as 


T h = C 



P AM Fs 

. Pm Fdh + Fs 


Mam 

Mmh 


(60) 
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In Eq. (60), if the term in square brackets is equal to 
zero,then T h will be constant. For the conditions given 
in Table 1 the magnitude of this term is equal to 0.2, 
which is quite small. It follows that whilst assuming lin¬ 
ear variations for the system pressure and masses of the 
expansion and compression volumes, the temperatures 
of the expansion and compression spaces may be 
assumed constant. 


7. Conclusions 

A numerical method is demonstrated for obtaining an 
approximate solution for the heat transfer in the regen¬ 
erative duct of free-piston Stirling engine with a recipro¬ 
cating internal wall. The method appears to be more 
suitable for free-piston type Stirling engines than other 
theories which assume constant gas and constant metal 
temperatures at a location in the duct with time. Another 
advantage of the present method is the inclusion of the 
effects of the gas that remains in the regenerative duct 
during the flow reversals. The assumption of a constant- 
pressure process in the system at an instant and constant 
heat transfer coefficient, based on the displacement, 
introduces errors in the prediction of the ineffectiveness 
of the regenerator. 


Appendix 

A simple method is developed to obtain a relation 
between the displacer and cylinder wall temperatures. 
The following assumptions are made to simplify this 
method: 

(a) There is no gas flow in the duct. 

(b) The longitudinal conduction is negligible in both 
the displacer and cylinder walls. 

(c) Heat transfer between the displacer and cylinder 
takes place in proportion to the temperature difference. 

The point Q on the cylinder wall with coordinate z is 
exposed to the temperature profile A B in Fig. AT The 
point R with coordinate z + dz is exposed to the tem¬ 
perature profile CD. For there to be no net heat 
transfer from a point on the cylinder wall to the dis¬ 
placer, the average value of the temperature profile to 
which a point is exposed must be the same as its own 
temperature. By considering the shaded areas in Fig. A1 
and their effect on the average height of the profiles to 
which the points Q and R exposed, the following equa¬ 
tion is readily obtained 


d Ty 
d z 


(Td\z+S/2) _ (^d)( r+ S/2) 


IS 


(Al) 



Fig. Al. Displacer-wall temperature distribution seen by 
points on the cylinder wall. 



Fig. A2. Cylinder-wall temperature distribution seen by points 
on the displacer wall. 


This equation is the condition for the quasi-static 
case, where the wall temperatures fluctuate periodically 
with small amplitude. 

Similarly, using Fig. A2 the following equation can be 
obtained from the consideration of the temperature profile 
on the cylinder to which a point on the displacer is 
exposed. 


d Tj 
dx 



(xxS/2) 


(Ty) 


(x—S/2) 


/S 


(A2) 


The initial temperature distributions of the displacer and 
cylinder walls are obtained using these relations derived 
between the cylinder and displacer wall temperatures. 
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